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In this paper we present in detail the numerical solution of the conformally invariant wave equa- 
tion on top of a fixed background space-time corresponding to two different cases: i) 1-1-1 Minkowski 
space-time in Cartesian coordinates and ii) Schwarzschild space-time. In both cases we use hyper- 
boloidal constant mean curvature slices and scri-fixing conformal compactification, and solve the 
wave equation on the conformal space-time. In the case of the Schwarzschild space-time we study 
the quasinormal mode oscillations and the late-time polynomial tail decay exponents corresponding 
^— V ■ to a mass-less scalar field. We also present general formulas to construct hyperboloidal constant 

^sj ' mean curvature slicings of spherically symmetric, static, space-times in spherical coordinates. 
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I. INTRODUCTION 



fj ' A common problem in numerical relativity is the global treatment of asymptotically flat space-times, specially 

qh| related to the implementation of boundary conditions on an artificial time-like boundary and the measurement of 

' ■ gravitational radiation by detectors located on time-like trajectories. Most of the initial value problems formulated 

{^Tjl are based on a Cauchy space-like foliation with the hypersurfaces approaching spatial infinity i". Then, boundary 

conditions must be supplied at an artificial boundary, for instance, absorbing boundary conditions [l| . Recent efforts 

point to the use of hyperboloidal foliations of the space-time, because asymptotically the spatial slices reach future 

null infinity ■J^'^ and for asymptotically flat space-times, the gravitational radiation is well dcflncd at such boundary 

[3, 3 J so that the measurement of the gravitational radiation can be observed by detectors near or at that boundary. 

[~^ I More common efforts are related to matching techniques of Cauchy slices with slices approaching future null infinity 

t^^ ■ [349, which has already been applied to the binary black hole problem Q. On the other hand, it is desirable to 

C^ I compactify the space-time foliated in this manner in order to work on a numerical domain that contains J^^ , for 

h — . which it is necessary to compactify the spatial coordinates in an appropriate way such that the space-time is regular 

f^ ' there, at future null infinity. At present time, foliations with boundaries at future null infinity are being used for 

^^ . various applications, like the solution of perturbation equations @t131 and the study of tails 11 



The idea behind scri-fixing conformal compactification is the following: i) use hyperboloidal slices that reach ^^ , ii) 
compactify the spatial coordinate, iii) rescale the metric with an appropriate conformal factor. We proceed by defining 
, , ■ a new time coordinate t = t — h{y), where {t,y) are the original time and space coordinates This transformation has 
rN ] the advantage that keeps the time direction invariant. That is, regardless of the choice of the h{y), the time-like 
S . killing vector has the same form in both coordinate systems. A compactifying coordinate is introduced in the form 
y = Q, where fl = fl{y) is a non-negative conformal factor that vanishes at ,y~^ [12|. Then the physical metric 
g becomes singular at the boundary of the compact coordinate. In order to remove the singularities the metric is 
rcscalcd g = fl^g with the conformal factor so that the conformal metric g is regular at the boundary J^"*". On the 
other hand, for the construction of h{y), we restrict in this paper to the case of constant mean curvature (CMC) slices 

In this manuscript we present the solution of the wave equation for the extended Minkowski and Schwarzschild 
space-times described with CMC slices, construct the foliations and compactification, and solve the wave equation 
on top of the resulting space-times with a twofold goal, on the one hand we verify as an academic exercise that the 
behavior of the wave function is consistent with the properties of the particular space-times constructed here, we are 
explicit in the construction of the numerical solution of the wave equation in one spatial dimension and illustrate the 
role played by the gauge. On the other hand, we study in particular the wave function on top of the Schwarzschild 
space-time, and track the quasinormal mode frequencies and tail decay rates of the amplitude of the wave function, 
which is also the case of a test mass-less scalar field evolving on the space-time. 

Once the space-time has been foliated with hyperboloidal slices, the spatial coordinate has been compactificd and 
the metric rescaled, it is necessary to solve a wave equation that is physically meaningful. Thus the wave equation to 
be solved has to be a conformally invariant one. Considering the conformally rescaled metric g = i^^g, the following 



identity applies for an n-dimcnsional space-time \\A 



n~2 , ^ 

4(n-l) 






= 0, (1) 



which means that the operator in brackets is invariant under conformal transforniations and is the one to be considered 
as the conformahy invariant D'Alanibertian operator. In this expression R and R are the Ricci scalar in the conformal 

and physical space-times respectively, = ^ is the conformal wave function, <j> the physical wave function and □ is 
defined by 

n0 = -La^[y^5'^-a,0], (2) 

where g^i, is the conformal metric and g its determinant. For the cases presented in this manuscript, the conformally 
invariant wave equation for the H-l Minkowski space-time is □</> = 0, whereas for the Schwarschild space-time (n ~ 4) 
the conformally invariant wave equation reads Ocfi — -rRcj) = 0. 

We solve these two equations as initial value problems in each case using a first order variable formulation. In the 
case of the solution of the wave equation on the Schwarzschild space-time, the equation also corresponds to a test 
mass-less scalar field, which is related to quasinormal mode oscillations of the black hole and study the polynomial 
tail decay of the amplitude of the scalar field; in such case we use mass-less scalar field and wave function indistinctly. 

This paper is organized as follows. In section II we solve the conformally invariant wave equation on top of the 
1-1-1 Minkowski space-time, and provide diagrams illustrating the results. In section III we solve the conformally 
invariant wave equation on top of a scri-fixing conformally compactified version of the Schwarzschild solution; we 
show the quasinormal mode oscillations and the tail decay of the scalar field, and compare with previous results. In 
section IV wc present general formulas useful to set up slicings with constant mean curvature and scri-fixing conformal 
compactification for spherically symmetric and static space-times in spherical coordinates; wc show as an example 
the 3+1 Minkowski space-time in spherical coordinates. Finally in section V wc summarize the main results. 

II. 1+1 MINKOWSKI SPACE-TIME 
A. Foliation and scri-fixing 

We choose the coordinates of the space-time x'^ = (i, x) and the line element 

dS^ = -dP +dS:'^, (3) 

where t, x e (— oo, +oo) x (— oo, +oo). In order to construct the hyperboloidal slices, we perform a time transformation 
t = t — h{x). Using this new time coordinate, the metric reduces to 

ds'^ = -dt^ - 2h'dtdS: + [1 - h'^]dS:^. (4) 

where h' = ^. In order to have a simple interpretation of the injection of hypersurfaces into the space-time we 
compare this expression with the usual ADM-like metric 



ds 



^2 



(-a^ + 72/32'! dt^ + 2(3^^dtdx + ^^dx^, (5) 



and we read off the lapse function a = . "'" , and then the future-pointing unit vector normal to the hypersurface 

\/l-h'^ 



n'^ is 



n^ = (-5,0) = ( ^===,0 1 , 



^r^^i 



2 



^'^^^^•TT^)' '«' 



with which it is possible to calculate the mean extrinsic curvature of the slice k — V^n'' ~ —^d^{^—gn^^): 

~k - di ( ,^' \ , ^ki= ( ^^' \ + C, (7) 

where the second equality is obtained after assuming k is constant, and finally C is an integration constant. This 
equation can be manipulated to give a complete description of the slices: 

kx 



h' = ^ => h{S:) = Va^+¥, (8) 

^l + (fcx)2 

where a = ~. The space-time in these new coordinates is described with hyperboloidal constant mean curvature slices. 
In order to compactify the space-time at the future null infinity boundary it suffices to define a new spatial compact 
coordinate x defined through x ~ f^. A convenient choice of the conformal factor is il — 1 ~ x^, because we want 
the two asymptotic ends x — ^ ±c» correspond to a; = ±1 in the compactified coordinate. This choice has been used 
for instance in J8| to construct the solution of the perturbation equation of charged wormholes. In this way, using 
X ~ x/n with fl — 1 — x^, h{x) = ^/a? + x^, implies h' ~ , ^ ^ , and finally the conformal metric reads 

ds^ = -il-xydt^-^J^^=£]^dtdx 



v/a2(l-a;2)2-)- 
a2(l + x2)2 



x^ 



a2(l -x2)2 +x 



dx'. (9) 



In order to find the gauge functions, this conformal metric has to be identified with the standard 1-t-l metric of type 
([5]), where we remove the tildes because we deal with the conformal metric. The identification implies the following 
gauge functions: 

rj.2 

2 ( /I 2n2 



^ _^ya2(l_a;2) + ; 



1 

„2n I 

2 



a-^(l — x^Y -\- x'^ 
which are the ones used to describe the space-time on top of which we solve the wave equation. 



B. The wave equation 

We solve the conformally invariant wave equation. For this case n = 2 and i? = 0, the wave operator is conformally 
invariant, thus we solve the equation ^ 

n0 = -La^[V^g^''a,0] = o 



for the metric ^ [TU|) as an initial value problem. 

It is straightforward to show that assuming a = a{x), /3 = /3(a;) and 7 — 7(2;), the wave equation can be written as 
a first order system of equations given by 

fa 

dti) = (9x — TT + /?V' 

V7 

V 7 



9x0 = ^, (11) 



where tt 



':df 



-(3dx(t> and ^A = dx<t>- This system of equations holds as long as 07 ^ 0. The third equation is 



a definition of "0 but it is also a constraint of the system that has to be satisfied. The value of the wave function 

— TT + P4' and can be integrated in time at the same time as 
so that the wave equation can be cast in a balance law form 

i9tu + A(9^u = —dx{A)vi, where 



is obtained from the definition of tt. that is dt(f) = 
TT and V'- We define the state vector u = (7r,-0)-^ 



P a/7 
a/7 /3 



(12) 



and the characteristic speeds of the system are given by 



dx 
'dt 



-/3± 



7 



(13) 



These eigenvalues are real and distinct and in fact correspond to a complete set of eigenvectors. Then, this system of 
equations is symmetric hyperbolic, and together with initial data for ip and tt it corresponds to a well posed initial 
value problem. We define such problem in the domain x £ [—1, 1], t € [0, 00). We choose initial data corresponding 
to an initially time-symmetric wave function with Gaussian profile. In particular we use cj) ~ e~^ /(oi) ^ -0 = Q^^ and 
7r = 0. 

We solve this problem using a finite differences approximation, with a method of lines using second order stencils 
along the spatial direction and a third order Runge-Kutta time integrator. The question about the boundary conditions 
7r(— l,t), '(/'(— 1,0: ^(1:^) and 7r(l,i) remains, and the answer is that there is no need to apply boundary conditions 
(e.g. radiative boundary conditions) because we notice that in the continuum limit at a:: = — 1 the characteristic speed 
corresponding to the mode moving to the right is zero (A+ = 0), whereas the characteristic speed of the mode moving 
to the left at a; = 1 is also zero (A_ ~ 0), and then no signals are propagated into the domain of the problem. 

The result is as expected, the initial Gaussian pulse will split into two pulses with half the initial amplitude and 
opposite directions, and the speed of propagation of each pulse depends on the coordinates chosen. The results of 
the evolution for three different values of the mean curvature k are shown in Fig. [TJ where the effects of the different 
gauge choices arc illustrated. In Fig. [5] we show the self-convergence of the solution, and the convergence of the 
constraint dx(\> — tp, which validates our numerical results. 
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FIG. 1: Evolution of the wave function for a — 0.5 (left), a = 1 (middle) and a = 5 (right) or equivalently k = 2,1,0.2 
respectively. These results are consistent with the space-time diagrams in Fig. 1121 where the light cones seem to get narrow 
near the boundaries for big values of a (small values of the curvature) which explains why the pulses when they approach the 
boundaries slow down and shrink. 



C. Space-time and conformal diagrams 

In order to better illustrate the structure of the space-time we construct the light cone structure and conformal 
diagrams, for which we first estimate the radial null geodesies by solving the condition ds = 0: 



dt 

dx 



x{l +a;2) 



± 



1 + x' 



(l-a;2)Va'(l-a;')'+a;' [l ~ x^) 
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FIG. 2: Convergence tests for the solution of the wave equation witli k = 2, 0.2 or equivalently a = 0.5, 5 on the left and 
right panels respectively. We use three resolutions Ax = 0.002, 0.001, 0.0005. (Top) We show the order of convergence Q of 
the constraint C = Tp — 8x4) to zero, defined by 1/2 (C) [using medium resolution]/ L2{C)[using high resolution] — 2*^ because 
the ratio between resolutions is 2. In the insets we show the L2 norm of the constraint for the three resolutions. (Bottom) 
Continuous lines represent the difference of (p for the two first resolutions and x —lines represent the difference of for the two 
last resolutions multiplied by 4; the fact that the curves lie on top of each other indicates the second order self-convergence of </!> 
for the snapshots presented. In the case a = 5 (small curvature) we can see that the errors are not as small as in the other case; 
this is due to the fact that the pulses are being squeezed when approaching the boundaries, and then the resolution used is not 
enough to resolve the pulses, and then the errors are bigger, however we maintain second the convergence of second order. 



whose solution is 



t^tn 



— y^a'^{l — x^)2 + x"^ ± X 



1 



(14) 



The result is shown in Fig. |3l What can be learned from this figure is that for small values of k the light cones become 
narrow near a; = ±1. The implication is that the pulses of the initial data propagating outwards will slow down and 
squeeze when they approach the boundaries in the case of small k. 

In order to obtain a better understanding of the nature of the slices, we construct the conformal diagrams of the three 
cases studied. The conformal diagram is constructed using the coordinates t = t{t,x) and x = x{t,x). Considering 
that h{x) 



\J a? + xP' and 51 



1 



the physical spatial coordinate reads x 



With this information one has 



the tilde coordinates given by: 




(1-^2) 



2~|2 ' 



l-a;2 



The results are depicted in Fig. 21 These diagrams are useful to interpret the effects of the gauge used. For the 
case of small k the hypersurfaces used for the evolution, even if they reach J''^ they are concentrated near spatial 






FIG. 3: Space-time diagrams for slicings with fc = 2, 1,0.2 or equivalently a = 0.5, 1,5. Continuous (dashed) lines represent 
null rays moving toward x = +1 {x — —1). These plots are constructed using (|14[1 with various values of to- 
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FIG. 4: Conformal diagrams with hyperboloidal slices for fc = 2, 1,0.2 or equivalently a — 0.5, 1,5. Continuous lines represent 
curves of constant t whereas dotted lines represent curves of constant x. 

infinity, whereas in the cases of bigger values of k slices are very separated near ^+ and allow the pulses arrive at 
the boundaries quickly. 

III. SCHWARZSCHILD SPACE-TIME 
A. Foliation and scri fixing 

We want to solve the wave equation on top of the Schwarzschild background in coordinates with two important 
properties: i) the wave is allowed to enter the event horizon of the black hole and ii) the wave is allowed to reach 
future null infinity. We then start with the Schwarzschild metric written as usual 



ds 



7.2 



2M 



dt 



n 



df 



?2 



2 A/ 



+ P(dr +sin^6ld(^^), 



(15) 



and proceed to construct a hyperboloidal foliation and a scri-fixing compactification. 

By introducing the change of coordinate t = t — h{f) in ([T5|) the line element takes the form 



ds^ 



1 



2M 
f 

2M 
f 



dt^ -~2h' 1 



2M 



r 
r 



dtdf 



df^ 



+ P(dr +sin^6idv3^), 

where h' — ^, and from which we read off the gauge quantities 



(16) 



7 
/3 = - 



f 



2M 



1 ^ I /i^. 

r 



(17) 



from which the unit normal vector to the spatial hypersurfaces pointing to the future is given by 



h' ( 2M\ 
7,^ 1-^ ,0,0 
7 V r J 



(18) 



Now, the extrinsic mean curvature at the initial time slice k ~ V^n^ reads 



k = ^df 



7 V f J 



(19) 



and can be integrated for constant k: 



where de derivative h' is 



h' 



kf^ 
~3~ 



C 



Ph' 



r) 



(20) 



(¥-g) 



fl-^^ 



(21) 



kr^ /^ 



(I - m) f4 



In this case, it is not easy to find h in a closed form as in the previous case in this paper, and in order to have a 
description of the slices one has to integrate this function numerically, however notice that h is not required but only 
its derivative. 

On the other hand, in order to perform the scri-fixing compactification we choose the conformal factor to be 
17 = 1 — r, and thus the compactifying coordinate is given by f = g. Then the Schwarzschild space-time using 
scri-fixing conformal compactification is described by the line element [13l | 



ds' 



where 



1- 



2Mn\ 



n^dt^ - 



r / 

„2 , „2/ ./i2 



P2(, 



P{r) = n^P{r) 



\ 



2{kr^/3-Cn^) 



dtdr 



dr^ +r^{d9-' + sm^edip^), 



Cfl^ 



2Mn\ 



fl?r 



(22) 



(23) 



The values of k and C are restricted in such a way that P(r) is real. Important features of these coordinates are 
that i) the event horizon is located at r = 2/3, ii) the slices penetrate the event horizon, iii) the slices do not avoid 
the singularity at r = f = 0, and iv) they reach the future null infinity. Condition (ii) allows the practice of excision 
inside the horizon, that is, a chunk of the domain is removed from the numerical domain which we know lies inside 
the horizon of the black hole, which means that the domain r < ri.r^c is removed and no need of boundary conditions 
at Texc are needed as long as the speeds of the characteristic fields are negative at that boundary. 



The gauge functions for this conformal metric are in this case 




(CI73 - ifca;3)2 



^(0 ^^^3 1 1^3 



^ ~ P{r) 



(24) 



B. The wave equation 

Following our convention, the conformal metric is given by ds^ = O^ds^, meaning that the physical metric is the 
one with the tildes. We decide to solve the wave equation using the conformal metric because of various reasons: 1) 
the slices are hyperboloidal and reach ^+ at infinity, 2) the spatial coordinate is compactified, and therefore the wave 
function reaches future null infinity at the boundary r = 1, 3) such boundary is null and there is no need to impose 
boundary conditions there, because the characteristic speed A_ = on that boundary. In order to use these benefits 
it is needed to solve the conformally invariant wave equation 

n-^R\cpT{t,r,e,^)=Q, (25) 

where R is the Ricci scalar of the conformal metric and □ = V^ V^ corresponds also to the conformal metric. In order 
to study non-radial modes, we separate (j)T{t,r,0,ip) = (j}{t^r)Yim{Q,^) with Yim the spherical harmonics. We solve 
numerically this equation considering a domain r G [vexci !]• We choose rf,xc such that it lies inside the event horizon 
and satisfies the need of P{r) being real. 

We solve ([25]) for ^{t,r,9,ip) as an initial value problem using a first order variable formulation. In terms of the 
line element given by the gauge (P^ the system of equations to be solved is 

dt^P = dr [ -TT + Pi: ] , (26) 



1 ^ /^„2/o_ , ",,^^ „.Y1d , ^(' + 1) 



,7 

dtTT = ^dr { r'{l3TT +-^)]^a-f[ ^R 
drCJ) = -0, 

where R — — 5- (r + m(2r — 1)). The system of equations is again symmetric hyperbolic and we solve it using initial 
data for tt and V' in the domain r S [rexc>l] and t G [0, cxd). As initial data we choose an approximately outgoing 
Gaussian pulse that on the initial hypersurface satisfies 

5t(r0)+A+9,(r0) = O, (27) 

where A+ is the outgoing characteristic speed given by \+ ~ — /? + — . The initial data reads as follow 

V'(0,r) = -2^1^^<t>{Q,r), (28) 

7r(0,r) = -^(0,r)-^i^fl-^ 



The numerical methods consist again in a finite differences approximation using the method of lines, however this 
time we use more accurate stencils in order to track small amplitudes related to the tail decay behavior of the wave 
function. We use now sixth order stencils along the spatial direction and a fourth order Runge-Kutta integrator in 
time. 



We choose an excision radius Vexc < 2/3 such that in the numerical domain P{r) is real for the values of C and k 
we have used. About the boundary at future null infinity r = 1, the characteristic speed of the mode moving to the 
left A_ — —f3 — af-f is zero, so that the wave propagates along the boundary and not toward the numerical domain. 
This implies that there is no need to apply boundary conditions there. 
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FIG. 5: Snapshots of the wave function for k = 0.4, 0.1, left and right panels respectively for the wave equation on the space- 
time given by p2|l . Again, slices with small curvature imply the outgoing pulse squeezes and slows down before it reaches 
J^"*", whereas for big curvature values the outgoing pulse reaches future null infinity quickly. The initial Gaussian is located at 
r = 0.8. The pulse splits into two pulses one moving toward scri plus and the other moving toward the horizon. We excise a 
part of the domain with re.xc ~ 0.6. 





FIG. 6: We show the light cones for the space-time (|22|l for the two values fc = 0.1 (left) k = 0.4 (right). It can be seen that 
for small k the light cones stretch near the boundary whereas for big values of k the cones are wide open at ,y^ . 



In Figiniwe show the effects of the gauge parameters chosen for the slices of the space-time for a spherical wave. In 
one case, when a large value of the curvature is used, the light cones near J^"*" are wide open and the outgoing pulse 
of scalar field arrives at such boundary without suffering a considerable deformation, whereas for smaller values of the 
curvature the pulse slows down and squeezes a little. To explain this effects we show in Fig. [5] the light cone structure 
of the space-time for two values of the extrinsic curvature of the slices. An important comment about the location of 
event horizons from Fig. [6] is also in turn. Event horizons are located in numerical relativity via tracking null surfaces, 
and the event horizon is calculated approximately as the surface from which the evolution of null 2-surfaces to the 
future diverge toward inside the horizon and toward J^"*" |15j . In this case, where hyperboloidal slices are compactified 
and the causal structure is maintained even after the conformal rescaling, we are sure that some of the null rays reach 
^+ located at r = 1; a fine tuned set of null rays in Fig. [6] would locate the exact threshold between those outgoing 
rays (rays in this case of spherical symmetry) really escaping to future null infinity and those that reach the excision 
boundary, being such threshold the location of the event horizon at r = 2/3. 
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C. Tails decay rates and quasinormal modes for the test scalar field on Schwarzschild space-time 



The study of wave tails in general relativity was first studied by Price and PuUin [ly, [iTf . Polynomial tail decay 
shows up after an oscillatory phase related to the quasinormal ring-down frequencies of black holes and can be regarded 
as being due to back-scattering of gravitational radiation off the curved background space-time. It is estimated that 
polynomial tail decay ^ t^ depends on the mode I that is being analyzed and on the trajectory of the observer, for 
instance p = —3 for time-like observers [1^ whereas for null observers at J^+ p = —2 [131 when I = 0, and in general 
p = —21 — 3 a.t finite distances and p = —I — 2 at ^+. 

Recent studies about the tail decay of a mass-less scalar field have been presented, which include the dependence 
on the location of the observer |18| . In this section we reproduce the decay rates and quasinormal frequencies for a 
mass- less scalar field with angular momentum on a fixed Schwarzschild space-time background for I = 0,1,2. In order 
to do so, we solve again the initial value problem for P6I27|28|) and study the amplitude of (j) for large values of time. 

D. Numerical results 

In order to study the oscillations and the tail decay it is enough to track the behavior in time of the amplitude of 
the physical scalar field (j), which reveals the oscillations and polynomial tail decay. Even though we are solving for 
the test scalar field on top of the conformal metric, the amplitude of the conformal scalar field cf) will show the same 
behavior in time as the physical scalar field, this is why showing the results for cj) suffices, since the conformal factor 
depends only on the spatial coordinate. We locate detectors of the scalar field at various positions in the domain, and 
relate the physical and conformal positions via the definition of the compact coordinate f = r/(l — r) and its inverse 
r = r/(l + f). 

The typical behavior of the wave function amplitude at late times is shown in Fig. [7] for various values of the gauge 
parameters. At initial times we have the quasinormal mode oscillations, after which the polynomial tail decay appears. 
We have also verified that the exponents are independent of a wide range of initial data profiles, both, amplitudes 
and widths of the Gaussian. 



C=2.0. K=OJ 
C=2.0. K=0.4 
C=5.0. K=1.0 
C=6.0. K=2.0 
C=9.0. K=3.0 




FIG. 7; We show the amplitude of the test scalar field in time, including the oscillating and the polynomial tail decay stages 
for different values of C and K measured at ^'^ on Schwarzschild space-time for I — 2 and Madm = 1. We found that 
each observer measures the same decay exponent p for different values of C and K and also the quasinormal mode frequencies 
measured are the same for different values of C and K. 



An important property that can be explored at this point is the frequency of the quasinormal mode oscillations. 
In Fig. [8] we show the amplitude of the scalar field for two different values of the ADM mass M , in one case we use 
M = 0.5 as usually done in previous papers, and in the other case we use the value M = 1 which is the one we use 
in all our further calculations. A fit of the oscillations shown in Fig. [8] with the function yle*'"+'^^-'* reveals that the 
frequencies for the mode / ~ 2, using M — 0.5 are lo = 0.749 and ujj = —0.178 which is consistent with previous 
studies [UllSl. 

In order to validate our numerical results we show in Fig. [HI an example of the convergence of the constraint tp ~ dr-cf) 
and snapshots of the self convergence tests of the scalar field. These tests indicate that the convergence of our results 
range between fourth and sixth order for quite a long time, which is consistent with the order of accuracy of our finite 
difference approximation and our MoL time-integrator. 
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M.„u= 1-0 





FIG. 8: We show the quasinormal modes and the tail decay rates of the wave function for different values of the mass of the 
black hole, measured an observer located at J^'^ for I = 2, K = 0.4 and C — 2. We found that the frequencies are different for 
differents ADM masses and the decay exponents are the same. The figure in the left panel is a zoom of the plot in the right 
panel for two cases that shows more clearly the dependence of the frequency on the mass of the black hole. It can be noticed 
also that the tail decay exponent is the same in all cases. 
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FIG. 9: We show the typical convergence tests for one of our runs, in this case for I — and resolutions Axi = 7.6 x 10^**, 
Aa;2 = Aa;i/2, Axs = Aa;i/4. In the left panel we show the order of convergence Q of the L2 norm of the constraint C = ^ — dr4> 
as defined in Fig 2, and in the inset the L2 norm of the violation of the constraint for the three resolutions. In the right panel 
we show the order of self-convergence Qsc of ({>, using the L2 norm of the differences between the value of the wave function for 
the various resolutions. 

We also measure the wave function for observers located at different positions from the event horizon. In Fig. [TUJwe 
show how the polynomial tail decay exponent depends on the location of the detector for 1 = 2, and similar behavior 
is found for the other modes. 

An efficient way of calculating the exponent p consists in using the fact that assuming (/)(x, t) ~ KV with K a 
constant, the exponent is given by p = ^°^^u\ as done in [l^, so that it is possible to track the time dependence of 
the exponent itself. In Fig. [Tl]we show the exponent p for I = 0, 1,2, measured at various positions and verify the 
bounds of the exponents to be between p = —21 — 3 for time-like observers and p = —I — 2 sX ^+. 

IV. GENERAL FORMULAS FOR THE SCRI-FIXING CONFORMAL COMPACTIFICATION 



We now develop the construction of the conformal compactification for a static spherically symmetric space-time 
described by the type of metric 
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FIG. 10: Amplitude of the wave function measured at different distances from the black hole's horizon for I — 2, which indicates 
that the tail decay exponent depends on the location of the detector. We use: A = 0.1, xo = 0.8 and a = 0.1. 



1=0 




1=2 



FIG. 11: We show the exponent of the polynomial tail decay for I = 0,1,2 measured by observers located at various distances 
from the black hole. From top to bottom the curves correspond to detectors located at f = oo, lOOOM, 250M, llOM, 50M, lOM 
and the remaining curves (for the case I — 0) correspond to detectors located even closer to the horizon. 



ds 



;;2 



a\f)dt^ + b^{r)df^ + f\d6^ + sin^ Odip^) 



(29) 



where h and a are assumed to be known metric functions, and proceed to construct a hyperboloidal foliation and a 
scri- fixing conformal compactification. By introducing the change of coordinate t = t — h{f) the Une clement takes 
the form 



ds^ = ~a^{r)dt^ - 2h'{P)a^{r)dtdr 

+ [52(f) - a^{f){h'{f)f] dr^ + r^{d0^ + sin^ edcj)'^), 



(30) 



where h' = ^. Comparison of this metric with the 3+1 metric ds"^ = {—a^ + ^'^j3'^)dt^ +2/3j^drdt + j^dr'^ + f^{d9^ 
sin^ 9dtp^) allows to read off the gauge and metric functions 



a{r) 



m 



a(r)b(r) 

7(f) 






(31) 



In terms of these functions, the unit normal vector to the spatial hypersurfaces pointing to the future is given by 



j{r) h' {f)a{f) 



a(f)b{r) ' ^(f)h{f) 



0,0 
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(32) 



Given a space-like slice, we can drag it along the time-like Killing vector. This will give a slicing where the time 
translation is along the Killing vector. Now, the mean curvature at the initial time slice, is given by 



y{f)b{ 



■zrdf 



f'^h'{f)a^{r) 
7(f) 



(33) 



and can be integrated for constant fc, which means the slices arc of constant mean curvature (CMC): 

k I r^a{r)b{r)dr - C ^ 
where now h'{f) is 



7(f) 



(34) 



h'if) 



[kl{f) - C]b{f) 



(35) 



a(f) J[fc/(f) - C]2 + a2(f)f4 



and 



/(f) 



f a{f)b{f)df. 



(36) 



In general, it is not easy to find /i in a closed form, so in order to have a description of the slices in general one has 
to integrate this function numerically. 

On the other hand, in order to perform the scri-fixing compactification, we define the compact coordinate r by 
f = h and we rescale the original metric using the conformal factor Q,. The space-time using scri-fixing conformal 
compactification is finally given by the conformal line element 

ds^ = -a^{r)^^{r)dt^ - 2h' {r)a^ {r)(n - rn')dtdr 



+ [b^{r) - a^{r){h'[r)f] [VL - r^'f 
+r^{de^+sm^ed^^), 



dr^ 



(37) 



where the functions a(r), b{r) and h'{r) arc functions of r. The conformal factor J7 determines various properties of 
the resulting conformal metric. 

For instance, for the Schwarzschild space-time, a^ = ^ = {1 — ^|^), with which our expression (|3S]), using il — 1 — r 
is reduced to the expression for h'{'r) obtained by Malec and Murchadha [Ij], and the final version of the conformally 
rescalcd metric is (P^ . 

A. The wave equation 

In general, from ([57)) we can read off again the gauge in terms of the ADM-like metric ds^ — {—o? -|- j^p^^dt^ -|- 
2fij^drdt -f ^'^dr^ + r^{d9^ + sin^ Odip^) and obtain the following gauge and metric functions 



fir) 



a^(r)h'{r){n-rn') 

7^(0 
{b\r)-a^{r)h''^{r)){n-rn') 



/\2 



(38) 
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which are the final metric functions used to solve the wave equation for a generic space-time. In order to solve (|25l) as 
a first order system, first order variables arc define as previously tt = idt(t> ~ -^Pdrij) and -0 — dr4'- Then the system 
to be solved is (pS)) with the gauge (|55)) for given initial data on a given spatial domain and appropriate boundary 
conditions. 

B. Example: 3+1 Minkow^ski space-time 

As an extra example, we present the case of the Minkowski space-time in spherical coordinates. We start with the 
line element for the Minkowski space-time in spherical coordinates x^ = (t, f, 9, 0) for the physical metric given by 

ds^ = -dp + df^ + r^{de^ + sin^ Odcj)^), (39) 

where t G (—00,00) and r € [0,oo). The integral (p6|) implies /(f) = f^/3, which substituted into (|35|) with C = 
implies 

Va + J' 

where a is a constant that is given in terms of the curvature by A; = -. It is also possible to integrate again for h and 
obtain a complete description of the slices with the height function h{x) = \/a?+x^, which completes the description 
of the conformal metric of the space-time. 

C. A particular case of the 3+1 Minkowski space-time 

In order to compactify this case we use a compact coordinate f = ^. Unlike the 1+1 Minkowski case, one only 
needs to regularize at r := 1 for which various choices arc evident: il = 1 — r or 17 = 1 — r^, in fact in |2l| for the 
solution of the Klein-Gordon equation with a 0"* potential, the conformal factor used is 57 = (1 — 7-^)/2, which is a 
conformal factor that shows a regular Ricci Scalar R at the origin. For the example developed here we choose the 
first choice and the conformal metric then reads: 

2r 

ds^ = -(1 - rfdP dfdT 

^a2(l -r)2 +7-2 



a2(l _r)2 +r 
from which we read off the gauge functions 



dr"^ + r^ide"^ + sin^ edcfP), (41) 



a' = {l-rf + i,. 
r 



r-2 



p = ^Va'(l-r)2+r2 

1" - 2n °'^2^ 2 - (42) 

a^(l — r)^ + r^ 
In order to visualize the space-time structure we calculate the null rays of the space-time which obey ds = 0, that 



IS 



dr (l-r)2 \^ ^a2(i_^)2 + 
whose solution reads 



dt 1 r \ 

± 1 , (43) 



^{a^ + l){T-lY + 2r-l ^ 1 
t — to — : ± 



1 



1 
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(44) 



The results are presented in Fig. [T2]for values a = 0.5, 1, 5 or equivalently A: = 6, 3, 0.6. The corresponding conformal 
diagrams are shown in Fig. I13[ where the injection of the slices in the space-time depending on the value of the 
extrinsic curvature is manifest. 



'^ ^^^ .^^''^^ ^^-"^ ^^-— — — — 





FIG. 12: Space-time diagrams described with hyperboloidal slices and fe = 6,3,0.6, or equivalently a = 0.5, 1,5. Continuous 
(dotted) lines indicate outgoing (in-going) null rays. As in the previous cases, the light cones close near the boundary for small 
values of the curvature. 
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FIG. 13: Conformal diagrams of the spherically symmetric Minkowski space-time in spherical coordinates, using hyperboloidal 
slices and conformal compactification with fi = 1 — r and various values of a = 0.5, 1, 5 or fc = 6, 3, 0.6. These diagrams are 
quite similar to those of the 1-1-1 Minkowski space-time, however notice that we are using a different conformal factor and 
that the curvature of the slices is different in terms of the parameter a. In fact, it has been found that a conformal factor 



equation 



)/2 is helpful to obtain a regular Ricci scalar at the origin, which is useful to solve the conformally invariant wave 



FINAL COMMENTS 



We presented the numerical solution of the wave equation on top of particular space-times using hyperbolidal slices 
and scri-fixing conformal compactification. We calculate the solution as an initial value problem using first order 
variables, and verified that the resulting constraints when the first order variables are defined converge, and also that 
the numerical solutions self-converge. 

We presented a detailed step- by-step construction of the foliation and conformal compactification for the l-fl 
Minkowski space-time. Following previous results wc constructed the case of Schwarzschild. We generalized in detail 
the scri-fixing conformal compactification for the case of a spherically symmetric static space-time, and provide a 
recipe to construct the adequate coordinates in that case. We illustrate each case with space-time and conformal 
diagrams in order to show the global structure of the space-time and how the slices inject into the space-time. 

The solution of the wave equation for the case of the Minkowski space-time in l-fl, is in fact the solution on the 
physical metric, because the wave operator is conformally invariant. In the Schwarzschild case -four dimensional- we 
want to stress that we solve the wave equation on top of the conformal space-time background. We studied the case 
of an outgoing pulse and tracked the late-time behavior of the amplitude of the wave function, and were able to study 
the quasinormal mode frequencies for various masses of the black hole and modes I = 0,1,2, and the exponents of the 
tail decay. We fitted the frequencies found previously and showed that the tail decay exponents obey the restrictions 
studied in the past. 
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